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A high statistics calculation, performed at P — 5.74, 6.00 and 6.26, enables us to study the variation of the 
leptonic decay constants fp of heavy pseudoscalar mesons with the lattice spacing a. We observe only a weak 
a dependence when the standard \/2k normalization is used for the quark fields, whereas application of the 
Kronfeld-Mackenzie normalization induces a stronger variation with a. Increasing the meson mass from l.lGeV 
to 2.3GeV this situation becomes even more pronounced. 



1. INTRODUCTION 

The prediction of the leptonic decay constants 
of the D and the B meson within the framework 
of lattice QCD is a challenging but also very deli- 
cate problem, since in the region of heavy mesons, 
the inverse of the respective masses comes close 
to currently reachable lattice resolutions. There- 
fore large discretization effects may contaminate 
the results . 

The question how to suppress these unphysi- 
cal contributions has been tackled from various 
sides Using meanfield arguments, Kronfeld 
and Mackenzie suggested that the replacement of 
the standard •\/2k normalization of Wilson quark 
fields by -^/l — 3k/4kc should signifcantly reduce 
the effects of finite lattice spacing. 

The present study takes a rather empirical ap- 
proach to the issue of finite a effects in the region 
of heavy mesons. We vary the lattice spacing in 
the currently accessible range and analyse the cor- 
responding variation of fp. Clearly a strong de- 
pendence of fp on a would indicate the presence 
of large discretisation errors, whereas a weak de- 
pendence would be associated with smaller con- 
taminations. In this sense our approach is per- 
fectly suited to judge on the efficiency of different 



quark field normalizations within a given a region. 

We will finally perform an extrapolation of fp 
to the continuum, assuming that its functional 
dependence on a is linear in the leading part, as 
suggested by the data. 

2. PREPARATION 

In order to visualize unambiguously the a de- 
pendence of fp, we have to take care that the 
finite a effects are not hidden in the statistical 
noise or distorted by incomplete groundstate pro- 
jection of the meson propagator and effects due 
to the finite size of the lattice. Therefore we have 
done our calculation with high statistics, keeping 
the errors of the raw data below 5%. We have 
varied the lattice size from about 0.7 fm to 2 fm 
and have smeared the quark fields with the well 
established Q Gauss like Wuppertal wavefunc- 
tion (n = 100, a = 4). In table |^ we display 
the lattice parameters together with the lattice 
spacing, taken from the stringtension a The 
influence of the finite lattice extension on fp was 
checked by comparing the results at different lat- 
tice sizes and fixed lattice constant. We find that 
finite size effects are small once the lattice exten- 
sion becomes as large as 1.4/m. 
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Table 1 

Lattice parameters 



(3 = 5.74, a-' 


= 1.118(9) 


P = 6.00, a-^ 


= 1.876(19) 


/3 = 6.26, a-' 


= 2.775(18) 


Ns, Nt 


no. con f 3. 


Ns, Nt 


no.confs. 


Ns, Nt 


no.confs. 


4,24 


404 


6,36 


227 






6,24 


131 






12,48 


103 


8,24 


175 


12,36 


204 


18,48 


76 


10,24 


213 


18,36 
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3. RESULTS 

3.1. Finite a effects 

In figj^ we show the Icptonic decay constant fp 
as a function of a both in the ^/2k normalization 
(open symbols) and in the Kronfeld-Mackenzie 
normalization (closed symbols). The light quark 
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Figure 1. fp as a. function of a. Points connected 
only by dashed lines were not used in the extrap- 
olation. 



mass has been extrapolated to the chiral limit 
and we have interpolated between the results at 
adjacent heavy quark masses (c.f. table ||) in or- 
der to keep the meson mass (in GeV) fixed when 
the a dependence of fp{Mp,a) is investigated. 
Due to their small statistical errors, we have used 
stringtension measurements Iq] to relate our data 



to a physical scale. The renormalization factor 
Za was taken from perturbation theory [^J^, with 
an effective coupling g'^ — Zg^/ < trP^^ >, re- 
commended in ref. (Pp,y =1x1 Wilson loop). 

It goes without saying that the a dependence 
of fp must be different in the two normaliza- 
tions. Very surprisingly, however, fig.|l| shows 
clearly that - in contrast to the •\/2k normalized 
results - the variation with a becomes stronger 
and stronger with increasing meson mass when 
the KroMac normalization is used. This means 
that - at least in the displayed a and Mp range 
- the KroMac normalization does a bad job: In- 
stead of suppressing finite a effects it enhances 
them. 

In order to connect our results to the (physical) 
continuum, we followed the behavior suggested by 
the data in both normalizations and extrapolated 
linearly||to a = 0. As can be seen from fig.|l|,we 
obtain nice agreement of the results, although the 
KroMac normalization has induced considerably 
larger errors. 

3.2. Heavy mass extrapolation 

The most 'natural' scale for fp is /^r, since the 
uncertainty originating from the renormalization 
constant Za cancels out in this case. Lattice mea- 
surements of /tt are generally affected with large 
statistical errors and therefore we have decided 
to convert our results to this scale only after hav- 
ing performed the a —^ extrapolation of fp. 
To achieve this we have decoupled the extrapo- 



lations according to ^(a 0) 



/,/v^(a-*0) 



1 - 0.0248g 



1 — 0.1333(;^ for standard normalization and Z a 
^ in the case of KroMac normalization. 



^ Since the a dependence cannot be exactly linear for both 
normalizations at the same time, we have excluded those 
points from the fit where .^1 — 3 -\/2k > 1.6. 
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Table 2 



Decay constant and meson mass in lattice units. The light quark has been extrapolated to Kc- 



/3 = 5.74 




P = 6.00 




/? = 6.26 








Jp/Za 


Mp 




fp/ZA 


Alp 


Kh 


fp/ZA 


Mp 


0.06 


0.1197(102) 


2.502(13) 


0.10 


0.0873(17) 


1.498(10) 


0.09 


0.0437(32) 


1.579(19) 


0.09 


0.1629(52) 


1.871(13) 


0.115 


0.0983(18) 


1.197(7) 


0.10 


0.0486(31) 


1.375(14) 


0.125 


0.1890(33) 


1.205(7) 


0.125 


0.1038(18) 


0.995(5) 


0.120 


0.0609(33) 


0.965(10) 


0.140 


0.1907(38) 


0.904(6) 


0.135 


0.1085(19) 


0.780(7) 


0.135 


0.0689(29) 


0.636(6) 


0.150 


0.1829(53) 


0.684(5) 


0.145 


0.1032(31) 


0.551(2) 


0.145 


0.0711(24) 


0.382(4) 














0.1492 


0.0580(50) 


0.245(4) 



Although the 0(g^)-uncertainty in Za does not 
cancel out exactly if one first extrapolates and 
then takes the ratio, its effect should be roughly 
the same in numerator and denominator. To 
obtain the denominator of this ratio we used 
both our own data and the results quoted in 
refs. Since the a dependence of ji^l^fo 

is weak, a linear extrapolation to a = is well 
justified and leads to ^(a = 0) = 0.269(12). 

In figure | we display our final results^ at a = 
in the formF fp{l/Mp), together with our static 
value from ref. [H. The new data appears to de- 



points carry error bars of order 25% and therefore 
do not exclude a stronger variation in Mp. Given 
this situation we draw an error band that links the 
conventional results with the static point. The 
Mp dependence of the error band was chosen ac- 
cording to the ansatz /p = cq -t- + jp- ■ 
At the location of the B and D meson the er- 
ror band corresponds to the bounds 155Mey < 
Jb < 242MeV , IdOMeV < fa < 200Mey . 
It is evident from figure ^ that these bounds are 
strongly affected by the size and uncertainty of 
/stat ■ More work is necessary to obtain an accu- 
rate prediction for /b- 
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Figure 2. / as a function of 1/Mp. 



pend only weakly on Mp. Because of the var- 
ious extra- and interpolations however, the data 

^The extrapolation has been performed on the data in the 
norm since it involves smaller statistical errors than 
using the KroMac norm. 
V"=/PvWx(H^)6/33. 
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